Practality of the non-equilibrium stationary states of open 
volume-preserving systems: II. Galton boards 
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Galton boards are models of deterministic diffusion in a uniform external held, akin to driven peri- 
odic Lorentz gases, here considered in the absence of dissipation mechanism. Assuming a cylindrical 
geometry with axis along the direction of the external field, the two-dimensional board becomes 
a model for one-dimensional mass transport along the direction of the external field. This is a 
purely diffusive process which admits fractal non-equilibrium stationary states under flux boundary 
conditions. Analytical results are obtained for the statistics of multi-baker maps modeling such 
a non-uniform diffusion process. A correspondence is established between the local phase-space 
statistics and their macroscopic counter-parts. The fractality of the invariant state is shown to be 
responsible for the positiveness of the entropy production rate. 

PACS numbers: 05.45.-a,05.70.Ln,05.60.-k 



I. INTRODUCTION 

Studying the statistical properties of simple mechani- 
cal models with strongly chaotic dynamics helps under- 
standing the connection between deterministic motion at 
the microscopic scale and transport processes which oc- 
cur at the macroscopic scales. This is of particular im- 
portance with regards to the irreversibility of thermo- 
dynamics and specifically the dynamical origins of the 
positiveness of entropy production. 

Such a mechanical device was originally introduced by 
Sir Francis Galton in the form of an apparatus which pro- 
vides a mechanical illustration of the Gaussian spreading 
of independent random events The Galton board, 
also known as quincunx or bean machine Q, consists 
of an upright board with a periodic array of pegs upon 
which a charge of small shots is released. The particles 
are let to collide on the way downward, thus displaying 
a seemingly erratic motion through the successive rows 
of pegs, until they reach the bottom of the board, where 
they are stopped. 

Provided the actual dynamics are sufficiently chaotic 
and dissipative, one can idealize individual paths as 
Bernoulli trials, whereby every collision event results into 
the pellets hopping down to the right or left of the pegs 
with equal probabilities. The number of steps in the tri- 
als is then specified by the number of the rows of pegs 
in the board. Under such conditions, the heaps of shots 
that form at the bottom of the board are expected to be 
distributed according to a binomial distribution and thus 
approximate a normal distribution. 

Though Galton's board was intended precisely as a me- 
chanical illustration of this idealized model, the dynamics 
of the board are necessarily more intricate, in particular 
with regards to inelasticity of the collisions between pegs 
and pellets and the friction exerted by the board's sur- 
face on the pellets. However if the collisions between the 



pellets and pegs were perfectly elastic and the board fric- 
tionless, the energy of every individual pellet would be 
conserved along its path. As a consequence, the kinetic 
energy would increase linearly with the distance separat- 
ing the pellet position from the top of the board, where 
one can assume it was released with a specified velocity, 
which, for the sake of specializing the motion to a fixed 
energy shell, we assume to be equal in magnitudes for all 
pellets. Such a conservative Galton board is also referred 
to as idealized. 

The remarkable property of conservative Galton 
boards is that a pellet's motion is recurrent, which is 
contrary to what had until recently seemed to be a 
widespread consensus. In other words, however far a pel- 
let goes in the direction of the external field, and con- 
sequently however large its kinetic energy becomes, it 
will come back to the top of the board with probability 
one. This property was proved by Chernov and Dolgo- 
pyat [2(J, , who also showed, in accordance to previous 
heuristic arguments and numerical studies, that the pres- 
ence of the external field affects the scaling law of posi- 
tions and velocities so that a pellet's speed scales accord- 
ing to v(t) ~ t 1 !' 1 and its coordinate x(t) ~ t 2 / 3 . They 
further found exact limit distributions for the rescaled 
velocity i _1 / 3 u(i) and position t~ 2 / 3 x{t). 

Galton boards and related models have attracted much 
attention in the statistical physics community. In par- 
ticular, Lorentz gases, which describe the motion of in- 
dependent classical point particles in an array of fixed 
scattering disks, have been the subject of intensive inves- 
tigations as models of diffusive transport of light tracer 
particles among heavier ones [!, 0, [E H, 0] • Lorentz gases 
have enjoyed a privileged status in the development of 
non-equilibrium statistical mechanics, which stems from 
the simplicity of its dynamics. By neglecting the recoil 
of heavy particles upon collision with the light tracer 
particles, one obtains a low-dimensional model that is 
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amenable to a proper thermodynamical treatment while 
it retains important characteristics of genuine many- 
particle systems. This model has been studied with 
mathematical rigor and, in particular, the existence of 
a well-defined diffusion coefficient has been proved rigor- 
ously under certain conditions [8|. Furthermore, in the 
last decades, and in the context of molecular dynam- 
ics simulations of non-equilibrium systems 0, Hof . sev- 
eral versions of the Lorentz gas model have been consid- 
ered, including the Gaussian thermostated Lorentz gas 
in the presence of a uniform external field [TTI ] , for which 
the Einstein relation between the coefficients of electrical 
conductivity and diffusion has been proved [T^. 

The reason for the initial success of the Lorentz gas was 
its use by Lorentz [§}, elaboratin g on Drude's theory of 
electrical and thermal conduction [13L [lij , for the sake of 
deriving the Wiedemann-Franz law, which predicts the 
temperature dependence of the ratio between heat and 
electrical conductivities in metals. In this framework, the 
computation of the electrical conductivity assumes that 
the external field is weak enough that the tracer parti- 
cle velocity magnitude is constant. Thus the diffusion 
coefficient is homogeneous and essentially given by the 
product of the particle's mean free path and (thermal) 
velocity. 

In a conservative diffusive system acted upon by an 
external field, the situation is different in that the exter- 
nal field causes the acceleration of particles and induces 
a velocity-dependent diffusion coefficient. Nevertheless 
such a system bears strong analogies with the field free 
diffusive case. 

It is our purpose to investigate this analogy by com- 
paring the statistical properties of Galton boards to that 
of periodic Lorentz gases. The latter were studied in a 
first paper where we discussed the fractality of the 
non-equilibrium stationary states of open Lorentz gases 
under flux boundary conditions, i. e. a slab of finite 
extension with its two boundaries in contact with parti- 
cle reservoirs with differing injection rates. Under such 
boundary conditions, the Lorentz gas sustains a steady 
current of mass which induces a constant rate of entropy 
production. 

In [l5| , we established the connection between this pro- 
duction of entropy and the fractality of the stationary 
states of open Lorentz gases. In this follow-up paper, we 
extend these results to Galton boards and related models. 
In particular, we develop a discrete random walk model 
that mimics the collision dynamics of Galton boards and 
associate to it a multi-baker map with energy, similar 
to models introduced in [l6|, [TtJ ■ The specificity of our 
model is that the transition rates vary with the sites' in- 
dices, reflecting the property of the conservative Galton 
board that deflection of tracers by the external field is 
more likely to occur when their kinetic energies are small 
in comparison to their potential energies. We derive 
the analytic expression of the non-equilibrium stationary 
states of this multi-baker map and show that its cumula- 
tive measures are characterized by nowhere differentiable 



continuous functions similar to the Takagi function of the 
non-equilibrium stationary state of the usual multi-baker 
map This allows us to compute the entropy associ- 
ated to such non-equilibrium stationary states and thus 
obtain an analytic derivation of the rate of entropy pro- 
duction, which, within our formalism, finds its origin in 
the fractality of the non-equilibrium stationary state, in 
agreement with the results presented in [f5j for the field 
free case. 

The paper is organized as follows. Galton boards 
are presented in Sec. HU The connection to the phe- 
nomenology of diffusion in an external field, described in 
Sec. Ill A[ is established for both closed and open systems, 
whose statistical properties are considered in Sec. IIIBI 
In Sec. Ill CI we discuss the occurrence of elliptic islands 
in the Galton board's dynamics, i.e. the stabilization 
of periodic orbits, and provide conditions under which 
we can assume the system to be fully hyperbolic. This 
regime is studied numerically, first under equilibrium set- 
ting in Sec. Ill Dl and then under non-equilibrium setting 
in Sec. Ill El where we demonstrate the fractality of the 
invariant measure. In Sec. IIIII we introduce the forced 
multi-baker map, which mimics the collision dynamics of 
the Galton board and analyze its statistics in Sec. IIII Al 
The entropy production rate of the non-equilibrium sta- 
tionary state is computed in Sec. IIII Bl We end with 
conclusions in Sec. HVl 



II. GALTON BOARD 

The Galton board is similar to a periodic Lorentz gas in 
a uniform field. We consider a two-dimensional cylinder 
of length L = Nl and height with disks E) n , < n < 
2N, of radii a, v / 3/4 < a /I < 1/2, placed on a hexagonal 
lattice structure. The centers of the disks take positions 



(nl/2,0), nodd, 
(nZ/2, ±V3Z/2) , n even. 



(1) 



where we identify the disks y = ±v / 3Z/2. The cylindrical 
region around disk B n is defined as 

l n = {(x,y) | {n- 1/2)1/2 <x < (n + 1/2)1/2}. (2) 

Thus the interior of the cylinder, where particles propa- 
gate freely is made up of the union U^' = _ Ar I„ \ HJ n . 

The associated phase space, defined on a constant en- 
ergy shell, is C = U^L_j V C 7l , where C„ = 8 1 ® [I n \ B„] 
and the unit circle S 1 represents all possible velocity di- 
rections. Particles are reflected with elastic collision rules 
on the border dC, except at the external borders, corre- 
sponding to x — 0, L, where they get absorbed. Points 
in phase-space are denoted by T — (x, y, v x , v y ), and tra- 
jectories by <E>'r, with $* the flow associated to the dy- 
namics of the Galton board. 

The collision map takes the point Y = (x,y,v x ,v y ) € 
dC to $T = (x',y',v' x , v' y ) S dC, where r is the time 
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that separates the two successive collisions with the bor- 
der of the Lorentz channel dC, and (v' x ,v' y ) is obtained 
from (v x ,v y ) first by propagation under the uniform ac- 
celaration until the instant of collision, and then apply- 
ing the usual rules of specular collisions. Given that the 
energy E is fixed, the collision map operates on a two- 
dimensional surface, which, when the collision takes place 
on disk n, is conveniently parameterized by the Birkhoff 
coordinates (<f> n ->£n), where 4> n specifies a generalized an- 
gle variable along the border of disk n, to be determined 
in Sec. Ill Dl and £ n is the sinus of the angle that the par- 
ticle velocity makes with respect to the outgoing normal 
to the disk after the collision. 

The external field is uniform and directed along the 
positive x direction, so that particles accelerate as they 
move along the axis of the channel, in the direction of 
the external field. There is no dissipative mechanism and 
energy is conserved along the Galton board trajectories. 

In this system, as opposed to typical billiards, the en- 
ergy, denoted E, can be both kinetic and potential. As 
the particle moves along the direction of the channel axis, 
it looses potential energy and gains kinetic energy, ac- 
cording to the energy conservation E — (v x + u?)/2 — ex, 
where e denotes the amplitude of the external field. Con- 
versely, the particle looses kinetic energy and gains po- 
tential energy as it moves in the direction opposite to the 
external field. 

Assuming E > 0, the boundaries of the system are 
placed at x — and x — L, reflecting the impossibil- 
ity for a trajectory to gain potential energy beyond the 
zero kinetic energy level. When E = 0, trajectories turn 
around at x > when the x component of the velocity an- 
nihilates, whereas when E > 0, depending on the choice 
of boundary conditions, particles can be either reflected 
or absorbed when they reach x = 0. 

The trajectory between two successive elastic collisions 
with the disks is now parabolic, according to x{t) = 
x(0) + v x (0)t + et 2 /2, whereas the vertical motion is uni- 
form y(t) = y(0) + v y (0)t. The amplitude of the external 
field e can be set to unity by an appropriate rescaling of 
the momenta and time variable: v — > v/\/e and t — > ts/e. 
Correspondingly, the energy has the units of length. 

We can thus write the velocity amplitude as a function 
of the x coordinate, 



v(x) = y/2(E + X ). (3) 

In particular, the velocity amplitude at x = is v(0) = 
\2E. We will assume that the energy takes half integer 
values of the cell widths I, so that the kinetic energy 
takes half integer values at the horizontal positions of the 
disks along the channel, i.e at x's which are half integer 
multiples of I. 

The system is shown in Fig. Q] with absorbing bound- 
ary conditions at x — and L. Note that trajectories are 
seen to bend along the field only so long as the velocity 
is small enough that the action of the field is notice- 
able. Otherwise the trajectory looks much like that of 



the Lorentz channel in the absence of external field. The 
time scales arc however different. 



A. Phenomenology 

One often reads in the literature that Galton boards, or 
equivalently periodic Lorentz gases in a uniform external 
field, do not have a stationary state. This is however a 
confusing statement since the existence of the stationary 
state has nothing to do with the presence of the external 
field. Rather, it is a matter of boundary conditions. 

Just as with the usual Lorentz gas, when an exter- 
nal forcing is turned on, a stationary state is reached so 
long as one specifies the boundary conditions. The rea- 
son for much of the confusion associated to this problem 
is, according to our understanding, that one cannot con- 
sider periodic boundary conditions along the direction 
of the field since they would violate the conservation of 
energy. One can however consider both reflecting and 
absorbing boundary conditions for the extended system. 
The nature of the stationary state, whether equilibrium 
or non-equilibrium, depends on the choice of boundary 
conditions. 

A phenomenological diffusion equation can be obtained 
for the motion along the axis of the cylindrical channel, 
which corresponds to the direction of the external field. 

In the presence of an external field, the diffusion pro- 
cess is a priori biased, so that the Fokker-Planck equation 
of diffusion reads 

d t v(x,t) = d x p(x)dxV(x,t) + M(x)r(x,t)}. (4) 

Here X denotes a macroscopic position, associated to the 
projection along the axis direction of a given phase-space 
region I„ of the Galton board, taken in the continuum 
limit. 

According to Einstein's argument, the diffusion coeffi- 
cient T>(X) is connected to the mobility coefficient Ai(X) 
by the condition that Eq. ((J) admits the equilibrium 
state V eq (X) as a solution which annihilates the mean 
current: 

V(X)d x V C(1 (X) + M (X)V cq (X) = . (5) 

At the microscopic level, letting T denote a phase point 
in 2d dimensions with velocity amplitude v and position 
x with respect to the direction of the external field, the 
equilibrium state is the microcanonical state i. e. 

Pcq(r) (X S(E - y + x\ , (6) 

Integrating this equilibrium phase-space density over 
cells C„ and taking the continuum limit I — > and 
n — > oo with the macroscopic position variable X — nl/2 
fixed, we obtain the macroscopic equilibrium density 
V eq (X), 

r cq (x)dx = iim / dF Poq (r), 




FIG. 1: Cylindrical Galton board with a non-vanishing external field and the energy E = 0. A trajectory is released at zero 
velocity at x — 0, and falls along the external field until it collides a first time with a disk. It then wanders around, coming 
back close to x — once, after which it moves further along the channel until it reaches the border at x = L. To compute the 
successive collision events, the numerical integration scheme uses an exact quartic equation solver based on the Galois formula 
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Identifying the length increments dX = I, and carrying 
out the velocity integration, we arrive to the expression 
of the equilibrium density 



V eq (X) =M[2(E + Xf 



(8) 



where Af is a normalization factor. Inserting this ex- 
pression into Eq. we obtain the relation between the 
mobility and diffusion coefficients, 



M{X) 



d-2 V(X) 
2~ E + X 



(9) 



The diffusion coefficient, on the other hand, is pro- 
portional to the magnitude of the position-dependent 
velocity, V(X) = \J2{E + X). This is a transposition 
of the corresponding result for the usual field free peri- 
odic Lorentz gas, where the tracer's velocity has con- 
stant magnitude. In the Galton board, given an en- 
ergy E identical for all the tracer particles, the veloci- 
ties V"(0) = \J2E at X — are identical for all particles, 
growing with X > 0, due to the uniform force of unit 
amplitude acting along that direction. We can therefore 
write 



P(W = P„\ 1+ J 



(10) 



Notice the normalization so chosen that the diffusion co- 
efficient at X — reduces to T>q. Equation (jTOj) can be 
thought of as a transposition of the argument by Machta 
and Zwanzig [191 ] who provided an analytical expression 
of the diffusion coefficient for the periodic Lorentz gas, 
based upon a random walk approximation. This approxi- 
mation indeed carries over to the Galton board. Provided 
energy is conserved, the velocity of a tracer particle in- 
creases as it moves along the direction of the external 



field. Thus, provided the periodic cells have sizes small 
enough that velocities remain approximately constant 
within each cell, the Machta-Zwanzig argument tells us 
that the diffusion coefficient is simply multiplied by a 
factor which accounts for the position-dependent veloc- 
ity. Hence the expression (jTOJ). 

Plugging Eq. (JTUJ) into ©, we obtain the expression 
of the mobility, 



M(X) = 



d-2 



-Do 



Ejl 



x_ 

E 



(11) 



Remarkably, the mobility coefficient vanishes for a two- 
dimensional billiard. In this case, the Fokker-Planck 
equation J2J therefore simplifies to 



d t V(X,t) = d x [V(X)d x V(X,t)] 



(12) 



An equivalent equation was derived by Chernov and 
Dolgopyat in [2(J. This is a diffusive equation without 
a drift and describes the recurrent motion of the two- 
dimensional Galton board trajectory at the macroscopic 
scale. In contrast, we notice that the Fokker-Planck 
equation Q associated to a three-dimensional version of 
the conservative Galton board has a non-vanishing mo- 
bility coefficient (JTTJ) and therefore retains a drift term. 

In the sequel we will assume E > so as to avoid the 
singularities that come with zero velocity trajectories. 

We notice, on the one hand, that reflection at 
the boundaries (RBC) induces an equilibrium state of 
Eq. I|12p with constant density, 

V{X) = 1, (RBC) . (13) 

Flux boundary conditions (FBC), on the other hand, viz. 



V(0): 

V{L) 



(14) 



admit the stationary state 

V(X) = C Q + C 1 ^/1 + X/E, (FBC). (15) 
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The coefficients Co and C\ are determined by the bound- 
ary conditions at X — and X = L, Eq. fjj 



Co 
Ci 



P-VT+ L/E-P+ 
y/l + LfE- 1 



y/T+L/E- 1 
In terms of V± , we can rewrite Eq. (fT5j> as 



= V- + {V+ - V-) 



(16) 
(17) 

(18) 



The corresponding local rate of entropy production is 
given according to the usual formula, by the product of 
the mass current (Q~9|) and the associated thermodynamic 
force H, 



diS(X) = [d x V{X)f 



dt 



T(X) 



(20) 



Given rates V- ^ V+, the current associated to the 
non-equilibrium stationary state is constant and, accord- 
ing to Fick's law, equal to 



J = -V{X)d x V{X), 
V V+ - V- 



2E y/l + L/E-l 



(19) 



Discretized Process 



The deterministic models we consider are to be an- 
alyzed in terms of return maps, which involves a dis- 
cretization of both time and length scales. We consider 
this problem in some detail, as it will be useful for the 
sake of defining a discrete process associated to the Gal- 
ton board. 



Let the discretized time and length scales be determined according to t — kr and X — nl. Collision rates are 
proportional to the velocity, which brings in a factor \Jl + nl/ E after we time discretize Eq. ([Tj 



TV E 



V(nl,kT + r) -V{nl,kr) 



(21) 



-2 \v(nl + l/2)V(nl + I, kr) + V(nl - l/2)V{nl - I, kr) ~ [D(nl + 1/2) + V(nl - l/2)]T(nl, fcr)} . 



We let 



A*n(fc) 



1 + %V{nl,kT) 
E 



(22) 



be the collision frequency on the Poincare surface at position X — nl, and introduce a diffusion coefficient associated 
to the discrete process, 



D(n) ee -Vinl) . 

It is convenient to set E = (2n$ + l)Z/2, for some positive integer no- 
Equation ([21) thus transposes to the evolution 



(23) 



Hn{k + 1) 



D(n+ 1/2) D(n- 1/2) 



In 
2n + l 



In 
2no + l 



D(n + 1/2) 
Vn(k) H (l n+ i{k) 



D(n- 1/2) 



1 + 



2(w+l) 
2n + l 



2(ra-l) 
2n + l 



Hn-i(k). (24) 



Written under the form 

H n {k+1) = s+_ 1 (i n - 1 (k)+s° n (i n (k)+s~ +1 iJ, n+1 (k) , (25) 

Eq. (|24| is seen to be the Frobenius-Perron equation of 
the Markov process 

f n - 1, s- , 

n — * < n, with probability , (26) 
I n+ 1, s+. 



As opposed to a symmetric random walk, the probabil- 
ities s~ , s° and s+ are asymmetric and depend on the 
site index, 





D(n-l/2) 




V /l +2n + l 


S+ 


£>(n+l/2) 






S° 

n 


= 1 - s n ~ 



(27) 
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In these expressions, n is assumed to be a positive integer, 
< n < N. From Eq. (fTU)) . the diffusion coefficient may 

be written D(n) = D 0y Jl^ ~ n /2; 

from which it follows that 



q-j-, where D a = t/1 2 V , 



8* = D Jl± 



2(n + n) + 1 



(28) 



It is straightforward to check that the stationary state 
of Eq. ([23)1 is independent of Do, and can be written 
under the form 



lim fi n (k) = 

k — KX> 



'2(n +n) + 1 



2nn 



Pn, (29) 



where P n is the discretized stationary state of the Fokker- 
Planck equation (TT2")) . 



Pn 



\/n + n + 1 + y/no + n 
= P n+ iVn + n + 1 + P„_iV«o + 



(30) 



We note that the latter equation implies that 
y/no + n(P n — P n -i) = a is constant. We can therefore 
write 

a 



Pn 



Pn-1 



\/n~o~- 



a > (n 

i=Q 



-1/2 



Pa 

Po+«{H* +no -m 



(31) 



where H 2 +n 



denotes the Harmonic number, Hn — 
E"=i J~ 1/2 & Letting n = N in Eq. ([nj), Nl — L, and 

writing the boundary conditions Pq = P_ and Pjy = -P + , 

i 

we obtain the expression of a, a = (P + — P-)/(H^ + — 



H 



no ). Therefore P n can be expressed as 

i i 
ff 2 — FT 2 



P n = P- + (P+ - P_) 



'n+n 



/T 2 — ff 2 



(32) 



The connection to the continuous case and, in partic- 
ular, to Eq. (I18|) is now straightforward. Indeed, the ra- 
tio of differences of Harmonic functions become integrals 
when / — > 



H 



n+n 



H 2 



ff 2 _ ff 2 



(33) 



y/E+l/2 y/E+31/2 



y/E+l/2+nl 



y/E+l/2 
r E+X 



I E ^ dx l 
Ie +L dx /\/x 



y/E+3l/2 ' ' ' y/E+l/2+Nl 
VE + X-\/E 

Ve + l-Ve ' 



where the limit assumes / — > with E constant and thus 
no ^> 1. In this case we have P n — > P(X = ni). 



C. Elliptic Islands 

Prior to turning to the stationary states of Galton 
boards, whether equilibrium or non-equilibrium, we men- 
tion the possible lack of ergodicity of the Galton board. 
The external field can indeed stabilize periodic orbits 
when the kinetic energy is not too large. Figure [5] shows 
such an example. In this case, elliptic islands co-exist 
with chaotic trajectories, as seen in Fig. [3] 




FIG. 2: (Color online) The external field induces a bifurcation 
such that the simple periodic orbit bouncing off two neighbor- 
ing disks at normal angles is replaced by two such orbits. In 
this situation where E = 0, which corresponds to a vanishing 
kinetic energy on the left border, one of these two periodic 
orbits is stable (to the left) and the other one unstable (to 
the right). The stability of the periodic orbit is quickly lost 
as E is increased. 

We notice that a mixed phase space is typically ex- 
pected in Hamiltonian chaotic systems -as is the case e.g. 
with the sine-circle map. This is an undesirable feature 
for our own sake. However the elliptic islands disappear 
if the energy value E is large enough. As it turns out 
of our numerical computations, E = 1/2 is already large 
enough. We will thus assume in the sequel that E large 
enough so the system is fully hyperbolic. 



D. Equilibrium Galton Board 

It is perhaps not widely appreciated that one can ob- 
tain an equilibrium state consistent with the presence of 
the external field. The reason for this is actually quite 
simple. Liouville's theorem implies the conservation of 
the volume measure, 



dT = dx dy dv x dv y , 
= v dE dtdipd^, 



(34) 



where <p and £ are defined to be the angle along the disk 
and sinus of the outgoing velocity angle measured with 
respect to the normal to the disk. 

We remark that because of the factor v that multiplies 
the volume measure in Eq. ([34]). the pair (0,£) are not 
canonical variables anymore. Indeed the position along 
the cylinder axis varies with the angle coordinate 4>, so 
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FIG. 3: Elliptic islands around the stable periodic orbit 
shown in Fig. [2] (Top) disk on the upper left corner; (Bot- 
tom) central disk. The phase-space coordinates used here (not 
the appropriate Birkhoff coordinates) are <j>, the angle around 
the corresponding disk, and £ the sinus of the velocity angle 
measures with respect to the normal to the disk. 



that the velocity v depends on <f>. The appropriate gener- 
alized angle variable conjugated to £ can be determined 
accordingly (36|. 

Introducing the index n, referring to the nth disk, 
whose center has position x = (n — 1)1/2 along the cylin- 
der axis, the velocity at angle <j> along disk n is 



v n {4>) = y/2[E+{n-l)l/2 + acos, 
= y (no + n)l + 2a cos 4>, 



(35) 



The canonical coordinate conjugated to £ is therefore ip n , 
such that d^ n = v n (4>) d(j>, 

I , ,n o Jo "n(0)# 
tj) n {(p) = 2tt-^ 



Jo* v n {<j))d4> 



2 2fT 

(n+no)^+2cr 



(36) 



E ( 2 (n+n )*+2 CT ) 

where E denotes the elliptic integral of the second kind, 
E(<j>,x) = J o 0/2 Vl-xsin 2 6d6, and E(at) = E(tt/2,.t) 
is the complete elliptic integral. As seen in Fig. 21 the 
difference between ij) n and 4> decreases rapidly as nl in- 



creases. Note that a is assumed to scale with I so that 
tp n does not actually depend on I. 
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FIG. 4: (Color online) Difference between the generalized an- 
gle coordinate ^>n(^>) and </>, here computed for a = 0.44/ , 
no = 0, and n = 1, . . . , 20. Larger differences occur at smaller 
n, where the effect of the external field is most noticeable. 



Let us consider a closed Galton board of length L = Nl 
(2N + 1 disks) , with reflecting boundaries at x = and 
x = L. This is an equilibrium system. More precisely 
the invariant density associated to each disk is uniform, 
as verified in Fig. [5] The distinctive feature however is 
that the time scale changes with the disk index n, r(n) ~ 
l/v n . Thus particles move faster with increasing n, but 
correspondingly they make more collisions so that their 
distribution is uniform in time. 

From the average count of collision events of disk n, we 
obtain the collision frequency, which, when multiplied by 
the local time scale (this amounts to dividing it by the 
velocity v n evaluated at the center of cell n) yields the 
average density P n ~ V(X n — (n — 1)1/2). This quantity, 
shown in Fig. [51 is indeed found to be almost constant, 
thus confirming our reasoning. 



E. Non-Equilibrium Galton Board 

A non-equilibrium stationary state of the Galton board 
can be achieved much in the same way as with the open 
Lorentz gas studied in [l5| , by assuming that a flux of tra- 
jectories is continuously flowing through the boundaries 
which are let in contact with stochastic particle reservoirs 
at x — and x = L. 



p(r,i) 



x=0X 



P±- 



(37) 



In analogy to the field free case, the invariant solution 
of the Liouville equation compatible with the boundary 
conditions (|37p is given, for almost every phase point T, 
by 



P(T) = P- 



P+ 



y/2(E + L) - V2E 



(38) 



8 



V* -! 




FIG. 5: Invariant density associated to a closed Galton board 
of size L = 10Z, with reflection at the boundaries x — and 
x — 101, and energy E — 1/2. This is an equilibrium system 
as reflected by the uniformity of the phase portraits. 
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FIG. 6: Equilibrium stationary density of the closed Galton 
board obtained for a channel of length L = 1, with 27V + 1 
disks, N — 25. The solid line is the constant equilibrium 
density V(X) = 1. 



v(T) 



'2E- 



-T(T) 



dtv(&r) 



Here p(T) is written in terms of the change in veloc- 
ity amplitude, given by v(T) = y/2[E + x{T)] at the 
corresponding horizontal position x(T). Thus v{T) = 
[x(T)]/v(T) and, provided the change in velocity between 
successive collisions is small, we can write JZ w(<i> t r)(it ~ 
[a;(<i> r r) — x(T)]/v(T). Hence, denoting by the time 
separating the (k — l)th and kth collisions and by tk — 
^2j=i T fe the time elapsed after k collisions, we have 



p(r) 



p- 



p- 



V / 2{E + L) - V2E 



v(T) - V2E 



K(r) 

E 

k=l 



a;($-**T) - x^-^T) 



(39) 



This approximation becomes exact when the number of 
cells in the system is let to infinity, in which case K(T), 
the number of collisions for the trajectory to reach the 
boundaries becomes infinite. Therefore the invariant 
state is 



p(r) 



p- 



p- 



E 

k=l 



_ v(T) 

^2{E + L) -V2E 

ar($-**r) - a;($-* fe - i r) 

u($-**r) 



>2E 



(40) 



so that the fluctuating part of the invariant density be- 
comes singular. This is analogous to the field free case 
discussed in [l5l ]. 

We compute this quantity numerically from the statis- 
tics of the Birkhoff map of the Galton board, using a 
cylindrical Galton board similar to that shown in Fig. [U 
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with external forcing of unit magnitude in the direc- 
tion of the cylinder axis, letting the particles have en- 
ergy E = 1/2. The particles are thus injected at x = 
with unit velocity at random angles and subsequently ab- 
sorbed upon their first passage to either x — or x = L. 

The computation of the collision frequency at disk 
n, averaged over the phase-space coordinates yields the 
quantity /i„, Eq. (|29p . which, after dividing by the mod- 
ulus of the velocity at that site, is converted to P n , the 
stationary solution of the Fokker-Planck equation (fT2l) . 
Here, we have 



T(X n 



c„} 



dTp(T) 



(T+ - V. 



^{e + l)-Ve 



.(41) 



The results of this computation are presented in Fig. [3 
and compared to Eqs. (fT5|) and (|3"2")l . The agreement with 
both discrete and continuous solutions is excellent. 
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FIG. 7: (Color online) Non-equilibrium stationary density of 
the Galton board obtained for a channel of length L = 1, 
with 51 disks (N — 25). Two solid lines are shown which are 
barely distinguishable, corresponding to Eqs. (|18p (Red) and 
(f32|) (Green) with E = 1/2 and thus n = 0. 



The histograms displayed in Fig.[H]show the fluctuating 
part of the invariant phase-space density computed in 
terms of the Birkhoff coordinates (ip n ,£,)i Eq. |36|) . The 
fractality of these graphs is much like that of the graphs 
of the open Lorentz gas, see [l5j]. The differences are 
indeed too tenuous to tell. As with the closed Galton 
board though, the distinctive feature is that the collision 
rates increase with the cell index with the amplitude of 
the velocity. 

To further analyze the fractality of the stationary state 
of the non-equilibrium Galton board and its relation to 
the phenomenological entropy production, Eq. ([20]) . we 
introduce in the next section an analytically tractable 
model, which generalizes the multi-baker map associated 
to a field-free symmetric diffusion process, so as to ac- 
count for the acceleration of tracer particles under the 
action of the external forcing. 



III. FORCED MULTI-BAKER MAP 

A time-reversible volume-preserving deterministic pro- 
cess can be associated to Eq. (|25|) in the form of a 
multi-baker map with energy defined on the phase space 
(n, [0, l n ] x [0, ln]) ne z, where each unit cell has area Z 2 = 
a n l 2 , a n = ^[1 + (2n)/(2no + 1)], and the dynamics is 
defined according to 



-i In — 1 x In — 1 + 



iV 



B : (n,x,y) 



n, 
n 



x — s„ l r . 



< x < l n s n , 

In^n — ^ — 

1. ^ ^"1^" , T^K+ik + s° n+1 l n + s- +1 y)) , l n (l - 4) < * < In- 



(42) 



This map has two important properties. First, the areas 
of the unit cells are chosen to vary with the amplitude 
of the velocity, which ensures that the Jacobian of B, 

is unity. Second, B is time-reversal 



or 



an+is n+1 



symmetric under the operator S : (n, x, y) — > (n, l n — 
y, l n — x), i.e. S o B = B^ 1 o S, as is easily checked. 

Multi-baker maps with energy have been considered 
earlier [l6|, [T3|- The novelty here is to introduce n- 
dependent rates and s°, Eq. P5|) . We will assume 
Da = 1/2 in the sequel, so that, provided uq + n ^> 1/2, 



we can write 

± _ 1 1 1 1 1 

Sn ~ 2 4 2(n + n) + 1 ~ 16 [2(n + n) + l] 2 

Thus s„ = 1 — s+ — s~ is approximated by 

o _ 1 1 
S "~ 8[2(n + n) + l] 2 + -" 

which is vanishingly small. Therefore, when n + n is 
large, the dynamics of B reduces to that of the usual 



(43) 



(44) 
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FIG. 8: Non-equilibrium phase-space densities of the open Galton channel with a geometry similar to that shown in Fig. [T] 
with absorbing boundaries at x = and x — 1, and stochastic injection of particles at x — only. Disk 50 is the one before last. 
Black areas correspond to absorption at the nearby boundary. The color white is associated to injection from the left boundary. 
Thus hues of gray correspond to phase-space regions with mixtures of phase-space points which are mapped backward to the 
left and right borders. The corresponding overall densities are shown in Fig. [7] 
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FIG. 9: (Color online) Forced mutli-baker map (|42|) . The 
cells have areas a n l 2 , with coordinates (x,y), and labeled by 
a positive integer n (no = here) . The map divides each cell 
into three vertical rectangles. The left rectangle is mapped to 
the bottom horizontal rectangle in the left neighboring cell. 
Likewise the right rectangle is mapped to the top horizontal 
rectangle in the right neighboring cell. The middle vertical 
rectangle is mapped to the middle horizontal rectangle in the 
same cell, but are so tiny they are barely visible already for 
n = 1. 



multi-baker map, at the exception of the energy depen- 



dence which fixes the local time scales. 



A. Statistical ensembles 

An initial density of points Y — (n, x, y), p(r,0), 
evolves under repeated iterations of B according to 
the action of the Frobenius-Perron operator, which, 
since B preserves phase-space volumes, is simply given 
by p(T,k + 1) = p(B~ 1 Y ,k). In order to charac- 
terize the stationary density, p(Y) = limfe—^ p(Y, k), 
we consider the cumulative function p n (x,y,k) = 
f£ dx' Jo dy'p(n,x',y',k). Notice that p here refers to 
the statistics of the return map and therefore differs from 
the density p associated to the Galton board, Eq. (|40p . 
by a factor proportional to the local time scale. 

The identification of this function proceeds along the 
lines of Refs. [H, [H, |3Zj]. Under the assumption that 
the x dependence of the initial density is trivial, we can 
write p n (x,y,k) = x/l n p n (l n ,y,k). Letting < y < 
1, it is then easy to verify that fJ> n (l n ,y l n , k) obeys the 
functional equation 



Pn{ln,yln, k+l) = < 



s n+lA*n+l { l n+l, 7f *n+lj 

S," 



< y < s+, 

n +lPn+l{Ui+l, Ui+li k) + S^p n (l n , V g „ n l n , kj , S+ < y < 1 — S % 



(45) 



(!-*«)<!/< 1. 



In particular, letting y — l n , we recover 

p n (l n , l n , k+l) = s~ +1 p, n+1 (l n+ i,l n+1 , k) (46) 

+ s n/*n('n! 'in + S n- 1 Mn— 1 On- 1 ) ^rt— 1> ^) ) 

which is identical to Eq. (|25p with p n (k) = Pn(ln,lmk). 
Let p n denote the steady state of this equation, p n — 
lim/c^ 

oo Pn(ln, 'n, 

The steady state of Eq. (|4"5|) can be written under the 
form, < x, y < 1, 

p n (xl n ,yl n ) = lim p n (xl n , yl n , k) , 

= a;y/i n + 2a; (s~ +1 ^ n+ i - s+p n )F n (y) , 
= xy p n + ax F n (y) , (47) 

where we introduced the generalized Takagi functions F n , 
with a prefactor, a = 2(s~ +1 /^ n+1 — s+// n ), which, as in 
Eqs. (|50|) - (f52")) is easily seen to be independent of n: 



2(s n+1 p n+1 ~ s+p 71 



! 2(n + n_+l) 
2n + 1 



(P„+i-P„), 



+ 



P_ 



2n 



Z 



1 rrl/2 



1/2 ' 



P 



+ 



2>/E y/E + L-^/E" 



(48) 



where the limit holds when I — ► 0. Notice that the pref- 
actor is proportional to Z. Thus a is a small parameter 
in that limit. 

Substituting Eq. (|4"T)) into Eq. (|4"5j) . the generalized 
Takagi function F is found to satisfy the functional equa- 
tion 



F n (y) = (49) 

0<y<s+, 
s+ < y < 1 - s~, 
1-s" <y< 1. 



The boundary conditions are such that the density is uni- 
form at n = 0,N + 1, implying F (y) = F N+1 (y) = 0. 
Notice that this function reduces to the Takagi function 
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in the limit n, N — > oo, n <C AT. Indeed s M ,,s+ = 1/2, 
s oo = 0- Therefore Eq. (|4"T|) is similar to the correspond- 
ing expression obtained for the multi-baker map, see [HI]. 



y(x) + |T n+ i(3a;) , 

2' 

1 - y(x) + 3 T„_i(3a; • 



2). 



< x < 1/3, 
1/3 < y < 2/3, 
2/3 < J/ < 1. 



i. Generalized Takagi functions 

For the sake of plotting F n (y), it is convenient to con- 
sider the graph of F n (y) vs. y as parameterized by a real 
variable, < x < 1, defined so that 



Vn{x) 



(50) 



SnVn+l{?>x) , 

s+ + s°y n (3x 



1) 



+ s ° +s n yn-i(3s-2), 



and 



F n {x) 



WW _l_ s 



2-S 

l-t)n(a:) 



n+1 

S n Fn (3^ 

S+ 
n- 



F„ + i(3a;) , 
1), 

i(3x 



2). 



< x < 1/3, 
1/3 < x < 2/3, 
2/3 < x < 1, 



(51) 

< x < 1/3, 
1/3 < y < 2/3, 
2/3 < y < 1. 



sequences {loi, . . . , w^}, ujj e {0, 1, 2}, 1 < j < k. 

The graphs of F n (xk) vs. y n (xk) are displayed in 
Fig. [TO] for a chain of iV = 100 sites and k = 8, and com- 
pared to the corresponding graphs of the incomplete Tak- 
agi functions [181] . which can be obtained from Eq. (|5T]) 
by setting s+ = s~ = 1/2 and s° = 0, 



2/0*0 



|y(3x) . 



+ ±2/(3* -2), 



< x < 1/3, 
1/3 < x < 2/3, 
2/3 < x < 1, 



and 



TJx) 



(52) 



(53) 



In passing, we note that, on the one hand, Eq. (|52|l is 
a functional equation whose solution is the Cantor func- 
tion. On the other hand, the tri-adic representation of 
the incomplete Takagi functions, Eq. (|53|) . is many-to- 
one. Their graphs, T n (x) vs. y{x), are nevertheless iden- 
tical to those obtained using the usual representation of 
the incomplete Takagi functions. 

2. Symbolic dynamics 

By substituting the tri-adic expansion of x in Eqs. (|50p ~ 

(ED, x({cj ,...,uj k }) = E- = o^ 3 7 (j+1) > «>i e {0,1,2}, 

we obtain the following symbolic representations of 
points y in cell n, 



Vn({uo, ■ ■ • , w fc }) = 



(54) 



uj = 0, 
= 1, 



s i + «° +s n y n -i({ui, ■ ■ ■ ,Uk}) , = 2. 



Starting from 



The boundary conditions are taken so that y n (x) = 
yi(x), n < 1, and y n (x) = yN(x), n > N. As above, 
F n (x) = 0, n < 1 or n > N. 

Starting from the end points y(n, 0) = 0, y(n, 1) = 1 
and F n (0) — F n (l) = 0, 1 < n < N, we successively 
compute y n (xk) and F n (xk), 1 < n < N at points Xk = 
E)=i 3_iw i, where, for every fc > 1, there are 3 fe different we can write 



Vn ({w }) 



0, 



w = 0, 

= 1, 
w = 2, 



(55) 



y„({w , . . . ,Wfc}) = 

2/n({w }) + ({wi, . . . ,Wfe}) , 

= 2M({^o}) + s^ w °y„ + i_ Wo ({wi}) 

+4" li,0 4+ii Wo y«+2-"o-^i({ w 2, ■ ■ • , w fe }) 



E 

i=0 



i-1 

n 

3=0 



1 — 

5 n+j-w -...-^i-i 



(56) 



Substituting this symbolic dynamics into the expression 
of F n , Eq. ((49]), we write 



\Vn+i{{ui, ■ ■ • ,w fc }) + s n+1 F n+1 ({uJi, . . .,uj k }) , uj q = 0, 

F n ({wo,...,Wfc}) = { l/2 + «°F Tl ({wi 1 ... > a; fc }) ) wo = 1, (57) 

|(1 - ■ • ■ ,w fc })) + s+_ 1 F n ^ l ({uj l , . . .,u k }) , uj = 2. 



r 



Let Ay„(cjo, ■ • ■ ,Wfc) denote the height of a horizontal {wo, • ■ ■ , Wfc}- We have 

J/„(w , . . . , uJk) 

y n ({uo, ■ ■ ■ , Uk + 1}) - 2/n({wo, • ■ ■ , w fc }) , 



cylinder set of the unit square, coded by the sequence ^ , cu ■) — (58) 
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112 1 

4 2 4 



Xk 




113 1 

4 2 4 

Xk 




Oilli 

4 2 4 

X k 




Ollli 

4 2 4 




Oiili 

4 2 4 
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4 2 4 
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FIG. 10: (Color online) Comparison between the graphs of F n (xk) vs. y n (xk) (Blue) and the corresponding incomplete Takagi 
functions T„(xk) (Red) . Each curve is computed at 3 10 + 1 different points x, uniformly spread between and 1. Only 2 10 + 1 
correspond to different points in the graphs of T„. 
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where the notation y n ({u>o, • • - ,u)k + 1}) is literal when- 
ever ujk ^ 2. Otherwise y n ({u>o, . . . , Wfe-i, 2 + 1}) = 
2/„({w , ■ • ■ , Wfc_i + 1,0}) and we set y n ({2, . . . , 2, 2 + 
1}) = 1. We have the following identities 

Ay„(0,wi,...,w fe ) = s+Ay„+i(wi, . . . , w fc ) , 
Ay n (l,wi, . . . , w fc ) = s°Ay n (wi, . . . , w fc ) , (59) 
Ay n (2,a;i,...,Wfe) = s~Ay n _i(o;i,...,a;fc). 

Therefore 

Ay n (w ,..-,w fc ) = n s «+i-wo---Wi-i' ( 60 ) 

t=0 

I 



which is nothing but the probability associated to the tra- 
jectory starting at position n and coded by the sequence 
{w , • ■ ■ ,<■*}. 

Likewise, the measure of the cylinder set 
Ay„(cj , . . . , u>k) is 

Afi n (uj , ...,(j k ) = fi n (y n ({uj , . . .,w fc _i + 1})) 

= /i„Aj/„(w , ■ ■ ■ , Wfc) + aAF„(w , . . . , oj k ) , (61) 
and we have the following set of identities for AF n : 



AF„(0,wi, . . . , u k ) = ±Ay n+1 (ui,...,u> k )+s n+1 AFn+l(u;i,...,uJk), 
AF n (l,u u ...,u k ) = s°AFn(o^ 1 ...,w fc ), (62) 
AF n (2,u>i,...,ujk) = -^Ay n -i(uj 1 ,...,cok)+s+_ 1 AF n ^ 1 (u)i,...,oj k ). 

That is, 

AF n (u , . .. ,cj fe ) = ^(1 - w )Ay n+ i_ Wo (a;i, . . . , cj fe ) + s^i Wo AF n+ i_ Wo (wi, . . . , uj k ) • (63) 

I 



Notice that it is possible to solve this system recursively, 
starting from 



AF n (w ) 



1/2, w = 0, 
, luq = 1 , 
-1/2, uj = 2. 



(64) 



We thus have a complete characterization of the non- 
equilibrium stationary state of B, Eq. (|42[) . associated to 
flux boundary conditions. 



B. Entropy and Entropy Production 

We proceed along the lines of [13, [HJ to obtain expres- 
sions of the entropies and entropy production rates asso- 
ciated to coarse grained sets such as defined in Eq. ([60)) . 
As described in [l5j, the idea is that, owing to the sin- 
gularity of the invariant density, the entropy should be 
defined with respect to a grid of phase space, or parti- 
tion, G = {dTj}, into small volume elements dTj, and 
a time-dependent state ii n (dTj,t). The entropy associ- 
ated to cell C n , coarse grained with respect that grid, is 
defined according to 



S t G (C n ) = -J2»n(dV j ,t) 



log 



dTi 



(65) 



This entropy changes in a time interval r according to 
A T 5*(Cn) = ^(C„)-C T (C„), (66) 



where, in the second line, the collection of partition ele- 
ments {dTj} was mapped to {Q T dTj}, which forms a par- 
tition <E> r G whose elements are typically stretched along 
the unstable foliations and folded along the stable folia- 
tions. 

Following [111, and in a way analogous to the phe- 
nomenological approach to entropy production [39l |. the 
rate of entropy change can be further decomposed into 
entropy flux and production terms according to 

A^(C„) = A^(C n ) + A^(C„) , (67) 

where the entropy flux is defined as the difference be- 
tween the entropy that enters cell C n and the entropy 
that exits that cell, 

AX(C„) = Sl^ dr . } (C n ) - Sl^ dr]] (^C n ) . (68) 

Collecting Eqs. (fS6")) - (|68p . the entropy production rate at 
C n measured with respect to the partition G, is identified 

as 



ATSt. 



^{* T drj}(^ 



(69) 



This formula is equally valid in the non-equilibrium sta- 
tionary state. 

Given a phase-space partition into the 3 fc cylinder 
sets coded by the sequences w_ k = {wo, ■ ■ ■ ,u>k-i}, Wi £ 
{0, 1, 2}, as described in Sec. IIII At the fc-entropy of the 
stationary state Eq. (|4"7]) relative to the volume measure 
of cell n is defined by 



<Sfc(C„) 



~ A Mn(Wfe) 



log 



Ay„(w fe ) 



1 



(70) 
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By summing over the first digit, it follows immediately 
from Eqs. (|4"5)) and ([SO)) that the fc-entropy verifies a re- 
cursion relation, 



S k (C n ) 



-S n+1 Pn+l fog ■ 

-s~ +1 5fe_i(C„. 



-l)- 

-l): 



S+-lPn-llog- 



s° n S k ^{C n ) 



+s+_ 1 5 fe _i(C„- 
with the 0-entropy given by 

So(C n ) = -/i„log/i„ 



(71) 



(72) 



and boundary conditions 



Sfc(Co) = -p-logp_ 
5 fe (Cjv+i) = -/°+l°g/°+ 



(73) 



The fc-entropy can be computed based on the above re- 
cursion relation. However, in order to obtain the de- 
pendence of the entropy on the resolution parameter k, 
it is more useful to consider the expansion of Eq. ([70]) 
in powers of p n . Let us denote by uj k the sequence 

{uj q , . . .,UJ k -l} 



S k (C T 



-fi n A v 



Ay n (u k ) 



log n Tt 



li n lAy n (uj k ) 



-/j„(log/^„ - 1) -a(log/i„ - 1) VAF„(w fe ) - ~ — V [^"(^fc)] 



0{af 



(74) 



The second term on the RHS of this equation vanishes, since 

^AF„(^) = 0. 



(75) 



As of the third term on the RHS of Eq. (f74l) . proportional to a 2 , we have, using Eqs. (|59|) and (|62|) . 



Aim . E^tf ■ 



1 1 

4s+ 4s.. 



-A^ +1 (fc-l) + ^A 2 (fc-l) + 



(^n+l) 2 

S + 
•5 n 



•7—1 



4s r t 4s. 
fc-i 

E E 

1=0 7/1 , — jTJi—l 



?J = 



J = l "tl+J — 1— m 



r 



(»t.i) s 



^-i(A-l) 



4s + 

v " n+'t — r/i — . . . — 7]i 



4s 



n+i— r/i — ...-r;,- 



(76) 
(77) 



Substituting the expressions of the probability transitions 
from Eqs. g3}-([44]), Eq. 1(77} is found to be 



+ O{n + n )- 4 . (78) 



A*(fc) = E 



A: 



[Afn(^ fc )] 2 

Ay„(w fc ) 
12fc 2 - 9fc 



32(n + n ) 2 

The first term on the RHS of this expression, which is 
the only term that survives in the continuum limit where 
n + no ^> 1, is responsible for the linear decay of the 
/c-entropy, 

a 2 

Sfc(C„) ~ -/i n (log/i„ - 1) - - — k . (79) 

2/^n 



Indeed if follows from Eq. (f6"5|) that the /Ventropy 
production rate is here 



ATS k (C n ) = -[5 fc (C n )-S fc+ i(C„)], (80) 

T 



Using Eqs. (TO]), (J2SJ) and (gSJ), it is readily checked that 
this expression yields the phenomenological entropy pro- 



duction rate, Eq. (gDJ), AT S k 



diS(X = nl)/dt. 
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IV. CONCLUSIONS 

In this paper, we have considered the influence of an 
external field on a class of time-reversible determinis- 
tic volume-preserving models of diffusive systems known 
as Galton boards or, equivalently, forced periodic two- 
dimensional Lorentz gases. 

Though the particles are accelerated as they move 
along the direction of the external field, in the absence of 
a dissipative mechanism, the motion is recurrent, which 
is to say that tracer particles keep coming back to the 
region of near zero velocity. In other words, particles do 
not drift in the direction of the external field. Rather, 
forced periodic Lorentz gases remain purely diffusive in 
two dimensions, albeit with a velocity-dependent diffu- 
sion coefficient. Consequently, the scaling laws relating 
time and displacement are different from that of a homo- 
geneously diffusive system. The macroscopic description 
through a Fokker Planck equation is however unchanged 
since the mobility coefficient vanishes identically in di- 
mension 2. 

It will be interesting to investigate the behavior of 
three-dimensional periodic Lorentz gases in a uniform ex- 
ternal field. As our analysis showed, the mobility does 
not vanish in dimension three, so that the Fokker-Planck 
equation retains a drift term. Being inversely propor- 
tional to the tracers' velocity amplitudes, this drift de- 
creases with increasing kinetic energy. A cross-over is 
thus expected between biased and diffusive motions. 

As far as their statistical properties are concerned, 
Galton boards are essentially identical to the field-free 
periodic two-dimensional Lorentz gases. A closed sys- 
tem with reflecting boundaries relaxes to an equilibrium 



state with a uniform invariant measure. This is to say 
that tracers spend equal amounts of time in all parts of 
the system. Open systems with absorbing boundaries 
yield non-equilibrium states. Given constant rates of 
tracer injection at the borders, the system reaches a non- 
equilibrium stationary state which is characterized by a 
fractal invariant measure. 

The fractality of the invariant measure associated to 
the non-equilibrium state of such a system was estab- 
lished analytically for a multi-baker map describing the 
motion of random walkers accelerated by a uniform ex- 
ternal field. The computation of the coarse grained en- 
tropies associated to arbitrarily refined partitions yields 
expressions which depart from their local equilibrium ex- 
pressions by a term which decreases linearly with the 
logarithm of the number of elements in the partition. 
This term is responsible for the positiveness of the en- 
tropy production rate, with a value consistent with the 
phenomenological expression of thermodynamics. 
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